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In the Standard Model, a group of heavy fermions, e.g. top quarks, can collectively strongly 
affect the Higgs field and create relatively long-lived bound states. If there exist new generations of 
fermions with masses beyond 1 TeV, strong binding of several of them can make them lighter than 
even a single heavy fermion. Using the mean field approximation we find multi-fermion states with 
masses M ~ 5v^/N « 1.2 TeV, with — 2, 3 . . . being the total number of heavy fermions 
bound together, and v = 246 GeV the Higgs VEV. The experimental search for multi-fermions within 
the range of energies 2 — 3 Tev would either discover them, or suggest absence of new Standard 
Model fermions with larger masses. Possible implications related to multi-top states and baryonic 
asymmetry of the Universe are discussed. 

PACS numbers: 14.65.Ha, 14.80.Bn, 12.39. Hg 
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The possibility that there may exist new generations of 
fermions, leptons and quarks, inspires search for quarks 
of the forth generation, see Ref. [l| and references 
therein. Since masses of known fermions cover a wide 
range of energies, one can/should anticipate that there 
exist super-heavy fermions, with masses m beyond 1 
Tev. We show that if they exist, their strong interac- 
tion with the Higgs field produces bound states for sev- 
eral heavy fermions (called multi- fermions below). The 
multi-fermion masses are relatively light, starting from 2- 
3 Tev for arbitrary large mass of heavy fermions. There 
exists therefore an opportunity for detecting super-heavy 
fermions via much lighter multi-fermions. This effect 
stems from the fact that fermion masses are propor- 
tional to the Higgs field. We will see that the Higgs field 
is strongly suppressed inside the multi-fermion which 
greatly reduces its mass. 

Implications related to strong interaction of fermions 
with bosons, in particular with the Higgs, which can 
modify its VEV have been discussed in Refs. 0, H, 
i, i, i, 0, i, 1 MM MMM- There were devel- 
oped the MIT-bag [^SLAC-bag[i| , and bag- version of 
Q. The Freiburg-Lee model 0] allows one to consider 
the MIT and SLAC models as its limiting cases. Our 
starting point, as well as the important soliton-type ob- 
ject, which appears in our analysis, is close to that dis- 
cussed in , as well as in a number of mentioned above 
works. An advantage of the present work is the ana- 
lytical solution, which we first derive for a large mass 
771 > V and then verify its applicability for a wide range 
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of fermion masses and numbers of heavy fermions using 
numerical analysis. These results provide clear physical 
picture of the problem, and give reliable estimates that 
may be used in the future for search of multi-fermions. 
Recently the idea that the quark can modify Higgs in 
its vicinity was tested for a single top and topponium 
13 , [iJl but only minor effects were found. Frogatt et al 
15| provided simple Hydrogen-atom-like estimates for a 
system of iV = 12 top quarks (6 top and 6 anti-top for 
the lowest 6*1/2 orbital), which suggested that its binding 
energy can be large. Unfortunately, our more accurate 
mean field calculations [lB| indicated that the binding is 
nonzero, though small, only for massless Higgs, while ac- 
count of the Higgs mass, which is limited by the current 
experimental bounds, makes this 12 top-antitop system 
unbound. 

Using the Standard Model, consider N heavy fermions 
that interact with the Higgs field. Take the conventional 
unitary gauge, in which the Higgs field $ is represented 
by the real field ^ 



V2U 



(1) 



Here v is the VEV, which is achieved when ^ = 1. The 
Lagrangian for the system of the Higgs and fermion fields 
^ and ip reads {h = c = 1) 



-1) 



- 7/5(77^9^- TTlOV^. (2) 
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The important feature of the problem is that either (i) 
the fermion mass m, or (ii) their number N or (iii) both 
are presumed large. In all the cases the impact on the 
Higgs field is strong, shifting its value away from the vac- 
uum VEV. We describe them here using the mean field 
approximation, though we are aware of effects due to 
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(weak and strong) gauge field forces as well as correc- 
tions beyond the mean Higgs field, which include many- 
body, recoil, and relativistic retardation/radiative effects, 
which would be addressed elsewhere [T^ . The magnitude 
of these corrections is different for cases (i)-(iii) because 
of their different nature and depends on various small 
parameters such as v/m, mn/m, and The case of 

large N > 6 (many top quarks) needs, in particular, an 
extensive discussion, as many levels are occupied subse- 
quently. 

Adopting this approach we replace the fields in Eq. ^ 
by stationary wave functions ^ for the Higgs and ^ for 
the single-particle wave functions of fermions. Searching 
for the spherically symmetrical solution, take f = ^(r) 
and assume that all fermions occupy the same shell with 
total and angular momenta j and /, which is described by 
the same large F{r) and small G{r) components of the 
Dirac spinor. Using Eq.Q one writes the Hamiltonian 
H of the system 



H ^ 



dr 



y(e'^ + i-l(e^-l)^)(4-'-^) (3) 

(^2{F' + ^F)G + m^{F^ ^ G^)) 

Here x = ±{j + 1/2) for I = j zt 1/2, and normalization 
J^{F^ + G'^)dr = 1 is taken. From Eq.® one derives 
the mean-field equations for ^, F, G 



(TV - l)m F2 - G^ 



r ' 2 ^ ^ ^ ' Attv^ 
{e-m£,)F = -G' + {K/r)G , 
(e-t-m^G = F' + {>c/r)F , 



(4) 

(5) 
(6) 



where the eigenvalue e is presumed positive [18]. 

It is interesting that the system considered reacts to 
the presence of the large fermion mass in such a way 
as to eradicate its influence on its physical parameters. 
This phenomenon employs three steps. Firstly, the Higgs 
develops a node on a sphere of radius ro, S,{ra) — 0, 
so that it is positive outside the sphere, taking at large 
distances the classical value ^ = 1, but is negative inside. 
Secondly, the fermions use this node of the Higgs as an 
opportunity to be localized in its vicinity, on the surface 
of the sphere, so that their density inside the sphere is low 
(hollow sphere). As a result the term in Eqs.(l5l),(l6|) is 
suppressed. Thirdly, the fermion wave function is tuned 
to satisfy F^{r) ~ G^{r), which eliminates the term ^ m 
from Eq.ljll). Thus, the described configuration of fields 
suppresses m everywhere in Eqs.Q-Q. 

To justify this physical picture analytically consider 
the large fermion mass, m 3> w. Assume that some 
smooth function ^(r), which has a node at tq, is given. 
Search for the solution of the Dirac equation in the form 



S{r) 



F{r) = A{r) exp(-5(r) ) 
G{r) = B{r) exp(-S'(r)) 
1 



m / ^(r')dr' 

'ro 



mCiro) {r - ro)' 



(7) 

(8) 

(9) 
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FIG. 1: Thick and dotted lines - large component F{r) of the 
Dirac spinor numerically and analytically, thin and double- 
dotted lines - same for small component G{r), slashed and 
slash-dotted lines - same for Higgs ^(r); numerical data - 
from solution of Eqs.llll-Q, analytical data for fermions from 
Eas.(ff))-(|9)). (|lip . analytical ^(r) is found from the minimiza- 
tion of Ti. in (I13p as explained in the text. 



The large mass m here favors localization of F{r) and 
G{r) in the vicinity of ro, which justifies the last identity 
in From ([S])-® one finds equations on A, B. To ver- 
ify that A and B are smooth functions of r in the vicinity 
of r = ro (which guarantees that fermions are localized 
near r = ro) we expand the coefficient functions >c/r and 
^(r) in (HI)-® in powers of r — ro, taking first the lowest- 
order approximation x/r w x/r^, ^(r) « C'(''o)(r — ro). 
Then one finds the eigenvalue 



e = ->^/ro = (j + l/2)/ro > , 



(10) 



which does not grow with increase of m, and is positive 
provided k < Q. The corresponding functions A{r), B{r) 
turn constants, A{r) = —B{r) = const. Normalizing 
them using Eqs.®-® we find 



A = - B = [m^'(ro)/47r]i/4 



(11) 



Using higher-order expansion for x/r and ^(r) in powers 
of r — ro we verified that Eqs. pop and (|lip remain valid 
for V. Thus, Eas.([7|)- PT|) give analytical solution of 
the Dirac equations ®,® when m > v. Fig. [1] shows 
good agreement of the found analytical wave functions 
with results of numerical calculations discussed in more 
detail below. 

Using the found eigenvalue (fTO|) of the Dirac equation 
we can simplify the Hamiltonian in Eq.®, where the 
term in the second line equals Ne — N\x\/ro. It is con- 
venient at this stage to scale distances by the Higgs mass, 
r X — mjj r, and present H in the form 

H = — n , 



n = in 



2 4^^ ' 



X dx H 

xo 



(12) 
(13) 



Here ^ — ^{x), ^' = d£,/dx, and xq = mn ro > is the 
node, £,{xo) = 0. The second term in Eq. (fT5|) . in which 

z = -N xm^/v^ > , (14) 
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FIG. 2: The etalon function h{z) presents the energy of 
a multi-fermion in Ea. (|12|l via the dimensionless quantities, 
H{z) and z from Eqs. (fT5)) . (fT4)) . 

reproduces the second terra from Eq. ([3|) , which equals A^e 
(as was mentioned). The integration in Eq. p^ neglects 
the contribution of distances x < xq- This approximation 
proves convenient and at the same time accurate since ^ 
is small and smooth for x < xq, while the factor x"^ in 
the integrand produces strong suppression compared to 
the outer region x > xq [l9j . 

The energy H of the system is presented in 
Eqs. p^ . P^ as a functional of ^{x). Finding its min- 
imum we derive the final answer to the multi-fermion 
problem in terms of simple, dimensionless quantities H 
and z. To find this minimum we minimized firstly Ti. over 
^{x) keeping the node xo fixed. This required solution of 
the boundary problem, + {2/x) + i'^)/2 = 0, 

^(xo) = 0, i{oo) = 1. The result found for H. was conse- 
quently minimized over xq. The outcome is the function 
Ti = Ti-iz), which we calculated numerically, depicting re- 
sult in Fig. [2]in terms of the etalon function h[z)^ which 
is related to ^{{z) via 

n{z) = 2(27rz)i/2(l + /^(z)) . (15) 

The first factor here equals the lowest term of the expan- 
sion oili{z) in powers oiz^/'^, Hiz) 2 {2TTzy/^ + 0{z) 
(with Xq [z/(27r)]i/2), indicating that h{0) = 0. It 
can be shown that h{z) reveals the asymptotic behavior 
h{z) - for z — > oo, but within a wide range of z the 
fimction h(z) remains small, h{z) <ti 1, see Figl2l From 
Eas. p^ - (|15p we derive the analytic expression for the 
total mass Ma of the system of N heavy fcrmions, which 
compose the multi-fermion 

Ma^2{2TTN\K\)^/^v[l + h{N\>i\m'^/v'^)] . (16) 

Here the term ^ h{z) remains small for a wide range 
of values of N,j and mn, see FigH] Consequently, the 
scale of multi-fermion masses is defined mainly by the 
Higgs vacuum expectation value v = 246 Gev, which 
gives Ma « 5.01 {N\x\Y/^v = 1.23{N\x\y/^ Tev. 

Consider numerical solution of the self-consistent mean 
field Eqs.lU])-®. Proper formulation of this problem in- 
cludes boundary conditions, which for fermions have the 



conventional form, whereas for the Higgs field they read 

e'(o) = , a^) = 1 . (17) 

The first one suppresses the singularity in the term 
(2/r)^' in Eq.(I3|) at r = 0. In linear equations this sin- 
gularity is harmless, it is eliminated by the conventional 
scaling of the function by ~ l/r factor. For nonlinear 
equations this, and others tricks do not work, the singu- 
larity persists making ^(0) singular. The first condition 
in p?)) removes this nuisance, allowing ^(0) to be finite, 
as it should be from the physical point of view. 

Our numerical calculations were performed taking 
TOH = 100 Gev for j = 1/2, l=^0,N = 2, 3, 6 and 12 
Fig. [1] shows wave functions for = 3 and m = 10 Tev. 
The good agreement of analytical and numerical results 
supports the fact that fermions in the multi-fermion oc- 
cupy mostly the surface of the sphere, where the Higgs 
field develops a node being positive outside the sphere 
and negative inside. Fig. [3] shows the calculated mass M 
of multi-fermions for different N and m. It is presented 
as a ratio M/AIa, where Ma is the analytical result (|16p . 
In the region of small m the multi-fermion mass M falls, 
and for sufficiently small m the bound state disappears. 
With increase of m the mass M shows a smooth max- 
imum, after which it decreases, revealing a tendency to 
converge to the analytical result (HH]). Fig. [3] shows that 
for a wide area of variation of m and N the masses of 
multi-fermions are close to the analytical predictions of 
Eq.dini). 

Remarkably, the multi-fermion masses found in 
Eq. (jl6p do not grow with an increase of the heavy fermion 
mass TO. For j = 1/2 Eq. p^ predicts masses 1.7 x K2 
Tev where 1.06 < if 2 < 1.17 and 2.1 x K3 Tev where 
1.07 < K3 < 1.21 for multi-fermions constructed from 
two and three heavy fermions respectively; the coeffi- 
cients K2 and K3 estimate corrections produced by the 
second term in the brackets in Eq. (fTB)l . when the Higgs 
mass spans the interval 100-300 Gev. Fig. [3] shows that 
if the fermion mass m is not very large, then the masses 
of multi-fermions are slightly larger, by up to 25%. 
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FIG. 3: Masses of multi-fermions M vs heavy fermion mass 
m, M is scaled by Ma from Ea. (|16|) . Solid, dashed, dashed- 
dotted, and dotted lines - A'' = 2, 3, 6, and 12 respectively, 
j = 1/2, l = 0,mH = 100 Gev. 
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For possible experimental applications it is interesting 
that for a small number of heavy fermions, N = 2,3, the 
masses of multi-fermions belong to the interval of energies 
of 2-3 Tev. This is in contrast to heavy fermions them- 
selves, whose masses may be larger. This fact provides 
an opportunity for hunting for super heavy fermions in- 
directly, via discovery of relatively light multi-fermions. 
A closely connected topic is related to bubbles made of 
a large number of top-quarks, or other heavy fermions 
(or bosons like Z,W). Note that the derivative dM/dN 
for the multi-fermion decreases with N. This means that 
the lifetime for the weak decay of a single fermion in- 
side the bubble increases fast with N because of decrease 
of the released energy M{N) — M{N — 1). This makes 
"magic" multi-fermions with complete or nearly complete 
ls,2p3/2, ■■• shells relatively long-lived states. Note that 
for a very large N the summation over occupied shells 
can be carried out analytically and may be expressed via 
the following redefinition of the parameter z in Eq. (|13p . 
z « (2/3)A^3/2^-i/2 j^y^2 ^here is the number of 
different types of fermions involved (e.g. 3 colors for 



quarks). Then the mass is given by Eqs. (fT2|) . (|15p . 

Another intriguing implication is related to the bary- 
onic asymmetry created during the Big Bang. Although 
in principle all Sakharov conditions are satisfied in the 
Standard Model, numerically they are extremely restric- 
tive. Nevertheless, as near the electroweak phase transi- 
tion the bubbles with different values of the Higgs field 
are created, the top (and possibly heavier) quarks get 
bound to the zero-Higgs surfaces, leading to long-lived 
states and large local deviations from thermal equilib- 
rium. This effect may substitute for the "bubble walls" 
(much discussed in literature when the 1st order phase 
transition was still an option). The CP violation and 
baryon number asymmetry in the decay of these states 
presents an appealing problem worth studying. 
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